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Abstract 


Six different convective-stratiform separation techniques, including a new 
technique that utilizes the ratio of vertical and terminal velocities, are compared 
and evaluated using two-dimensional numerical simulations of a tropical 
[Tropical Ocean Global Atmosphere Coupled Ocean-Atmosphere Response 
Experiment (TOGA COARE)] and midlatitude continental [Preliminary Regional 
Experiment for STORM-Central (PRESTORM)] squall line. The simulations are 
made using two different numerical advection schemes: 4th order and positive 
definite advection. Comparisons are made in terms of rainfall, cloud coverage, 
mass fluxes, apparent heating and moistening, mean hydrometeor profiles, 
CFADs (Contoured Frequency with Altitude Diagrams), microphysics, and latent 
heating retrieval. Overall, it was found that the different separation techniques 
produced results that qualitatively agreed. However, the quantitative differences 
were significant. Observational comparisons were unable to conclusively 
evaluate the performance of the techniques. Latent heating retrieval was shown 
to be sensative to the use of separation technique mainly due to the stratiform 
region for methods that found very little stratiform rain. The midl'atitude 
PRESTORM simulation was found to be nearly invariant with respect to 
advection type for most quantities while for TOGA COARE fourth order 
advection produced numerous shallow convective cores and positive definite 
advection fewer cells that were both broader and deeper penetrating above the 
freezing level. 
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Introduction 


Squall lines are an important class of precipitating convective systems, often 
referred to as mesoscale convective systems or MCSs, that have long been the 
subject of study by the meteorological community. They occur both over land 
and water, and both in the tropics and midlatitudes. They can contribute a 
significant amount of rainfall to a particular region, produce severe weather, and 
impact the large scale environment. Squall lines consist of a linear group of 
thunderstorms or cumulonimbi usually along the leading edge of the system in 
the direction of propagation and, in the mature stage, an associated anvil 
containing a broad area of light precipitation. The precipitation in squall lines 
can generally be categorized as being either convective or stratiform in nature. 
There are several reasons for classifying the precipitation and its region of origin 
within the squall line. First, the mechanisms by which precipitation is formed 
are decidedly different in the two regions. The convective region is also 
generally associated with heavier precipitation rates. In addition, the latent 
heating profiles associated with each region are distinctly different from one 
another with each region having its own characteristic shape. The same is true 
for the mass flux and divergence profiles. Thus, it can be useful to distinguish 
the two regions for a variety of purposes. As such, a number of techniques have 
been developed over the years that are designed to separate convective systems, 
including squall lines, into their convective and stratiform components. The 
purpose of this study is to compare several of those techniques within the 
context of a numerical cloud resolving model in light of their potential 
applications. 

An early study on tropical squall lines by Hamilton and Archbold (1945) 
described a typical squall line as a line of cumulonimbus clouds that preceded a 
broad downdraft region associated with a precipitating anvil that emanated from 
the cumulonimbus towers. In an early study on midlatitude squall lines, 
Newton (1950) constructed a time-height cross-section of a pre-frontal squall line 
showing a subsidence inversion behind the line using upper-air soundings from 
the Thunderstorm Project. Zipser (1977) looked into the kinematic structure of 
the squall line downdraft region using aircraft and sounding data gathered for a 
Caribbean squall line. He put forth a conceptual model wherein the downdraft 
region consisted of narrow near-saturated convective-scale downdrafts and a 



broader unsaturated mesoscale downdraft that originated from mid-level inflow 
and was driven by rain evaporation. These two sources of downdrafts led to the 
formation of a stable layer behind the leading edge convection. Brown (1979) 
used a two-dimensional (2D) hydrostatic model with parameterized cumulus 
convection and simple microphysics to confirm this theory of an evaporatively 
driven mesoscale downdraft in the anvil region. 

Houze (1977) had similar findings about the downdrafts in his study of a 
GATE squall line system to that of Zipser (1977). Using GATE radar data he 
computed the anvil percentage of the rainfall to be about 40%. He simply drew a 
boundary around the squall line region (i.e., the convective portion) and treated 
the rest as anvil. He noted the squall line had discrete cumulonimbus elements 
that formed the leading edge. New elements formed ahead of the system and 
merged with it. Later as they became older they blended into the stratiform anvil 
region. The precipitation falling from the anvil region was stratiform in nature 
as evidenced by a distinct radar bright band. First explained by Ryde (1946), the 
bright band is associated with broad areas of frozen light precipitation falling 
through the 0 C isotherm (Battan 1973; Atlas 1990). 

Leary and Houze (1979) used radar data from several GATE cases to 
estimate cooling rates in and below the anvil bright band due to melting and rain 
evaporation. They found the horizontally uniform precipitation of the anvil 
region to contain rainfall rates on the order of 1 to 10 mm/h compared to 10 to 
100 mm/h in the convective cores. Gamache and Houze (1982) obtained 
quantitative divergence and mass flux profiles for both the convective and anvil 
regions by compositing radar and wind data for a GATE squall line. The 
convective region contained low level convergence and upper level divergence 
v/hile the anvil portion had mid-level convergence and divergence at both 
upper and lower levels. Net mass flux profiles were positive at all levels for the 
squall region while the anvil region showed positive fluxes aloft but downward 
fluxes below the base of the anvil. They separated the convective and anvil 
portions based on subjective assessment of the echo intensities, reflectivity 
gradients, the presence of vertical cores, and steadiness of the patterns. They 
found the squall-line echo to have an average width of 20 km. Churchill and 
Houze (1984), using aircraft data collected from a tropical cloud cluster during 
Winter MONF.X, found much higher ice concentrations in the convective 



region. The convective region was also dominated by rimed particles or graupel 
as opposed to aggregates in the anvil region. 

An important aspect of squall line research involves their influence on 
the larger scale notably in the form of heating and moistening. Using a 
compositing technique, Reed and Recker (1971) estimated the average structure 
of 18 disturbances that crossed through a sounding network in the Marshall 
Islands. They obtained an estimate of the diabatic heating effects above 900 mb 
for a composite easterly wave disturbance. The difference between the trough 
and ridge positions showed positive heating effects throughout the troposphere 
with a maximum near 400 mb. Yanai et al (1973) used a similar approach to 
estimate the apparent heating and moistening for tropical disturbances also 
observed in the Marshall Islands area. Their time averaged heating profile 
matched that of Recker and Reed (1971) fairly well with a single maximum in 
heating near 475 mb. These early estimates represented heating for an entire 
convective system. Houze (1982) used an idealized cloud system to put forth the 
first comprehensive estimates of the heating profiles associated with the 
convective and anvil regions by summing estimates of the cloud-scale terms. 
The convective profiles showed warming dominated by condensation 
throughout the depth of the troposphere. The anvil profiles had condensational 
warming in the middle and upper troposphere and cooling due to evaporation 
and melting in the lower troposphere. Johnson and Young (1983) computed the 
heating and moistening rates for tropical anvils using Winter MONEX data by 
computing the large-scale terms using soundings from a ship array when the 
array was covered mainly by mesoscale anvil clouds. Their results compared 
well with those of Houze (1982). Their anvil moistening profile showed 
moistening in the lower troposphere and drying aloft. Johnson (1984) 
partitioned the apparent heat and moisture source profiles of Yanai et al. (1973) 
into cumulus and mesoscale components by estimating the fraction of rainfall 
associated with the mesoscale anvil and using the anvil profiles of Johnson and 
Young (1983). The cumulus profiles were obtained as the residual. The resulting 
heating profile for the cumulus region matched fairly well with that of Houze 
(1982). 


Tao and Soong (1986) used a three-dimensional (3D) numerical cloud 
model with warm rain microphysics to simulate a GATE rainband. They 



produced heating and moistening profiles for the entire convective system 
similar to those from observational studies. Tao and Simpson (1989) added ice 
microphysics and a convective-stratiform separation technique based on 
Churchill and Houze (1984) to produce the first model simulated profiles of 
heating that distinguished the convective and anvil regions. These convective 
and anvil heating profiles were similar to those obtained by Houze (1982) and 
Johnson (1984). The simulations were based on a GATE composite case. 

Recent studies (Simpson et al. 1988; Adler and Negri 1988; Tao et al. 1993a) 
have indicated that separation of convective and stratiform clouds is necessary 
for the successful retrieval of surface rain and latent heating profiles via remote 
sensing. Frank and McBride (1989) concluded that differences in the total heating 
profiles between AMEX and GATE cloud clusters were due to differences in the 
fraction of anvil rainfall. Using similar logic, Tao et al. (1993a) put forth an 
algorithm to retrieve the total mean latent heating profile for a system or region 
based on the rainfall, the percentage of rainfall that was stratiform and 
appropriate profiles of heating associated with the convective and stratiform 
regions. A higher stratiform percentage resulted in a more elevated level of 
maximum heating. The convective and stratiform heating profiles used in the 
algorithm could be obtained from both diagnostic and modeling studies of 
convective systems from various geographic regions. They indicated that the 
stratiform percentage should be to within 10% for an accurate retrieval. 
Alexander and Cotton (1998) recently devised a new mesoscale parameterization 
scheme to accompany a modified traditional cumulus parameterization scheme. 
They used the Tao et al. (1993b) convective- stratiform separation technique to 
select mesoscale profiles for their parameterization lrom cloud resolving model 
simulations of MCSs. 

Recent cloud resolving modeling studies (e.g., Tao et al. 1993b; Xu 1995; 
Caniaux et al. 1994) have quantified the amount of stratiform rain as well as the 
heating and moisture budgets associated with the convective and stratiform 
legions of mesoscale convective systems using different techniques to separate 
the convective and stratiform regions. Current convective-stratiform separation 
for TRMM ground validation radars follows the algorithm of Steiner et al. (1995). 
Indeed the exact demarcation between convective and stratiform is somewhat 
arbitrary, but such separation has many useful applications. As such, it is 
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important to compare and contrast the various separation techniques that are 
available to numerical cloud models to determine whether similar conclusions 
are achieved and the magnitudes of any quantitative differences. This paper will 
help to address the issue by comparing five different convective-stratiform 
separation techniques that are currently being used plus an additional new 
method based on the premise that the terminal velocity of precipitation particles 
is large relative to the vertical velocity in regions of stratiform precipitation (e.g., 
Houghton 1968; Steiner el al. 1995; Houze 1997). These techniques are applied to 
two different cases: a well organized mid-latitude case (PRESTORM) and a well 
organized tropical case (TOGA COARE). Section 2 describes the model, the two 
cases, the separation techniques and two different types of numerical advection 
that will be used to further test the sensitivity of the model with regard to 
convective and stratiform differentiation. Section 3 compares the results of the 
numerical experiments in terms of rainfall, cloud coverage, mass fluxes, 
apparent heating and moistening, mean hydrometeor profiles, CFAD (Yuter and 
Houze 1995b) analyses, microphysics, and latent heating retrieval. Finally the 
summary and conclusions are given in Section 4. 

2. Numerical Experiments 

The model used in this study is the two-dimensional (2D) version of the 
Goddard Cumulus Ensemble (GCE) model. The GCE modeled flow is assumed 
anelastically balanced, filtering out sound waves by neglecting the local variation 
of air density with time in the mass equation. The cloud microphysics include a 
parameterized Kessler-type two-category liquid water scheme (cloud water and 
rain), and parameterized Lin et al. (1983) or Rutledge and Hobbs (1984) three- 
category ice-phase scheme (cloud ice, snow and hail/graupel). Short-wave (solar) 
and long-wave (infrared) radiation parameterizations as well as a subgrid-scale 
turbulence (one-and-a-half order) scheme are also included in the model. Details 
of the model description can be found in Tao and Simpson (1993) and Simpson 
and Tao (1993). 

Two different cases, a well organized tropical oceanic squall system 
observed during TOGA COARE (Tropical Ocean Global Atmosphere Coupled 
Ocean-Atmosphere Response Experiment) and a well organized midlatitude 
continental squall system observed during PRESTORM (Preliminary Regional 



Experiment for STORM-Central), will be used to test the various convective- 
stratiform partitioning methods and advection schemes described in Sections 2.1 
and 2.2. The June 10-11, 1985 PRESTORM case has been well studied (e.g., 
Johnson and Hamilton 1988; Rutledge et al. 1988; Tao et al 1993b; Tao et al. 1996). 
The PRESTORM environment contained a fair amount of instability but is also 
relatively dry when compared to the tropical environment. The thermodynamic 
profile used to initialize the model was a single sounding taken at 2330 UTC 
from Pratt, KS which was ahead of the newly forming squall line. The sounding 
is quite unstable with a lifted index of -5.37 and a CAPE of 2300 J/kg. The 
PRESTORM simulations were made using a modified shear profile [see Tao et al. 
(1993b)] and Lin et al. (1983) microphysics which are more representative of a hail 
environment usually associated with such strong continental systems. Observed 
reflectivities for this case were characteristically over 50 dBZ (Rutledge et al. 
1988). A stretched vertical coordinate was used in the model with. 31 grid points. 
The grid resolution ranged from 240 m at the lowest level to 1250 m at the top. 
The horizontal grid consisted of 1024 total grid points; the central 872 of which 
had a fixed 1 km resolution. The outer grids were stretched. Both long- and 
short-wave radiative parameterizations were used [see Tao et al. (1996)] though 
the majority of the simulation occurred during nighttime hours. There were no 
surface fluxes used in the model, and the time step was 6 seconds. The 
convective system was initiated using a low level cold pool. A cooling rate of 
0.01 K/s was applied over a period of 10 minutes. 

The February 22, 1993 TOGA COARE squall line has also been the subject 
of numerous studies. This system was observed with Doppler equipped 
instrumented aircraft (Jorgensen et al. 1997), and has been numerically simulated 
(Trier et al. 1996; Trier et al 1997; Redelsperger et al 2000; Wang et al 1996, 2000). 
The sounding used to initialize the model for this case comes from LeMone et al. 
(1994). It is a composite sounding that consists of aircraft data below 6 km and a 
combination of the 1800 and 2400 UTC Honiara soundings above 6 km. The 
surface conditions are assumed to be: Psfc=1006 mb, Tsfc=28.66 C, and Qsfc=23.2 
g/kg (Redelsperger et al. 2000). The resulting CAPE and lifted index values are 
moderately unstable at 1776 J/kg and -3.2, respectively. A low level westerly jet 
of 12 m/s is present near 2 km. The observed squall line propagated eastward 
and was oriented perpendicular to this jet feature. The squall line was initially 
linear but later evolved into a bow shape with a mid-level vortex along the 
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northern edge. Since this was an oceanic case, surface fluxes were included using 
the TOGA COARE flux algorithm put forth by Fairall et al. (1996). Additional 
information on the implementation of this algorithm into the GCE model can be 
found in Wang et al. (1996). The use of the TOGA COARE flux algorithm 
necessitates the use of a very fine first grid level. As such, a similar type of 
vertical stretching to that in the PRESTORM setup was used, but with the first 
grid level at 40 m. This in turn required a small time step of 5 seconds. 
Resolution at the highest model level was 1150 m. The horizontal grid followed 
that for the PRESTORM case but with a non-stretched resolution of 750 m and 
slightly smaller stretching. A modified version of the Rutledge and Hobbs (1984) 
3-class ice scheme was used which is characteristic of graupel physics rather than 
hail. The shear above the low level jet at 2 km was reduced to allow for the 
development of a better stratiform region. Both solar and long-wave radiative 
fluxes were included as well. A low level cold pool similar to the PRESTORM 
case was used to start the convection. In addition, low level mesoscale lifting 
was applied in the TOGA COARE case. This lifting had a peak value of 3.4 cm/s 
near 1 km and was applied over the first 2 hours of the simulation. Table 1 
shows some of the characteristics of the large scale environments associated with 
these two cases. 

2.1 Numerical advection schemes 


A second- or fourth-order horizontal advection scheme has generally been used 
in the model. However, using second-order or higher-order accuracy advection 
schemes can introduce some difficulties because negative values arise in the 
solution (Soong and Ogura 1973). This effect can be especially important in cases 
where the solution of the advection is used as input to nonlinear equations 
describing microphysical phenomena or inert tracers which can eventually lead 
to instability of the whole system (Smolarkiewicz 1983). The use of upstream 
differencing or other low-order schemes (Soong and Ogura 1973) would not 
produce dispersive "ripples" but would suffer from excessive numerical 
diffusion. Smolarkiewicz (1983) has reduced the implicit diffusion by using a 
second upstream step where a specifically defined velocity field leads to a new 
form of a positive definite advection scheme with small implicit diffusion. The 
positive definite advection scheme involves iterations and needs more 
computational resources. This scheme has been improved to include 



multidimensional applications (Smolarkiewicz 1984) and a non-oscillatory 
option (Smolarkiewicz and Grabowski 1990). 

The GCE model has implemented this Multi-dimensional Positive 
Definite Advection Transport Algorithm (MPDATA). Two specific solutions are 
performed in this study. The first one has all scalar variables (potential 
temperature, water vapor, turbulent coefficient and all five hydrometeor classes) 
use forward time differencing and the MPDATA for advection. The dynamic 
variables, u, v and w, use a second-order accurate advection scheme and a 
leapfrog time integration (kinetic energy semi-conserving method). The second 
more traditional approach uses leapfrog time integration and a fourth-order 
space derivative scheme for scalar and dynamic variables. 

2.2 Convection-Stratiform Partitioning Methods 

(a) Churchill and Houze (19841 

Churchill and Houze (1984) designed a method to partition convective 
and stratiform regions using radar estimated rainfall rates. Their technique was 
based on an assumption from a previous work by Houze (1973) that convective 
cells have peak rainfall rates that are at least twice as large as the surrounding 
background. Consequently, they identified convective cores as points whose 
rainfall rates were at least twice that of the surrounding 400 km2. For each core, 
the surrounding 150 km2 was also deemed convective. Finally, any rainfall rate 
greater than 20 mm/h (~40 dBZ) was also made convective. Originally, 
Churchill and Houze (1984) applied their method to radar data at a 3 km altitude 
with a 4 km horizontal resolution. In the GCE model, this technique is applied 
to surface rainfall rates and normally on a higher resolution grid reducing the 
size of the background area (2 points in each direction). Also, the 20 mm/h 
convective threshold is allowed to vary between 10 and 25 mm/h depending on 
the time of the simulation. This method is termed the "C&H method" hereafter. 

(b) Tao and Simpson (1989) and Tap et al (1993bl 


In this convective-stratiform partitioning method, convective, stratiform 
and non-surface precipitating regions are identified using the information from 



surface rainrate first (C&H method). Two additional criteria are applied which 
have been included to identify regions where convection may be quite active 
aloft though there is little or no precipitation yet at the surface, such as areas 
associated with tilted updrafts and new cells initiated ahead of organized squall 
lines (Tao et al. 1993b). Non-surface precipitation regions are considered to be 
convective if cloud water exceeds a specific threshold (i.e., the maximum of 
either 0.5 g kg-1 or half of the maximum cloud water content at that specific 
simulation time), or if the updraft exceeds a specific threshold (i.e., 3 m s-1) below 
the melting level. The presence of this amount of cloud water is a good 
indication of a saturated area (100% Relative humidity). For the region classified 
as stratiform by the C&IT method, additional criteria are considered to identify 
active convective cells aloft. First the updraft velocity is checked above the 
melting level. If it exceeds the maximum between 5 m/s and half of the 
maximum updraft value at that simulation time, the column is considered 
convective. Second cloud water and cloud ice content above the melting level 
are also checked. If their combined value exceeds the maximum between 1.5 
g/kg or half of their maximum combined content at that simulation time, the 
column is also convective. 

These thresholds are obtained by examining the modeled cloud structures 
thoroughly. However, these threshold values may need to be adjusted for 
different cloud systems as well as for the various stages of their life cycle. 

Cloudy regions ahead of the gust front are also classified as convective in 
order to achieve coherent areas of convective and stratiform rainfall. This 
method is termed the "GCE method" and has been adapted by Chin (1994) and 
Alexander and Cotton (1998). 

(c) Xu (19951 

A third convective-stratiform partitioning method was applied to a cloud 
ensemble simulation by Xu (1995). This method's idea is very similar to the one 
developed by Churchill and Houze (1984), except that the horizontal gradient of 
absolute air vertical velocity below the melting level rather than rainfall is used. 
Model grid points that have a vertical velocity twice as large as the averaged 
value taken over the four surrounding grid points are identified as convective 



cell cores. For each core grid point, the grid point residing on either side is also 
made convective. Furthermore, any grid point with a rainrate in excess of 25 
mm/h or an absolute vertical velocity in excess of 3 m/s is designated as 
convective regardless of the above criteria. The remaining columns where the 
total liquid water path exceeds 0.2 kg/m2 are then defined as stratiform. The 
stratiform columns are then further checked for shallow convection. Columns 
that have cloud water paths above 0.4 kg/m2, rain water paths below 0.1 kg/m2 
and mean upward motion all below the melting layer are made convective. 
However, in this study, any areas with surface precipitation that were not found 
to be convective by the first three criteria (i.e., vertical velocity gradients, vertical 
velocity thresholds, and surface precipitation thresholds) were considered to be 
the stratiform region. In addition, the surface rainrate threshold was not fixed at 
25 mm/h but was allowed to vary between 10 and 25 mm/h depending on the 
simulation time. One other modification involves the 3 m/s vertical velocity 
threshold. In the GCE model it is taken to be the minimum between 3 m/s and 
one fourth of the maximum vertical velocity at that simulation time. Again, 
these modifications consider the resolution of cloud characteristics during its life 
time. This method is termed the "Xu method". 

(d) Caniaux et_al {1994} 

In this method, a characteristic width or fixed number of model grid 
points (e.g., 20-40) are centered on the maximum surface rainrate and designated 
as the convective region. Remaining grid points with surface precipitation 
comprise the stratiform region. This method is termed the "constant area or CA 
method". In this study, the convective region was fixed at a width of 20 km, so it 
is referred to as the "CA20 method". 

(e) Steiner et al 11995) 


This method is a texture algorithm applied to radar reflectivities at a level 
below the melting band and is an extension of the earlier C&H method. 
Reflectivities are compared against a background average taken over a radius of 
11 km. Points that exceed the background average by a certain threshold are 
made convective. The threshold varies as a function of the background average. 
In addition, any gridpoint that exceeds 40 dBZ is made convective. Finally, for 
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each of the gridpoints identified as convective, a surrounding area which 
depends upon the intensity of the core gridpoint is made convective. The 
remaining rainy areas are then stratiform. 

(f) Vt-W Method 


In this study, a new method is introduced. This method is based on the 
basic premise that the terminal velocity of precipitation particles is large relative 
to the vertical velocity in regions of stratiform precipitation (e.g., Houghton 1968; 
Steiner el al. 1995; Houze 1997). As such, the magnitude of the terminal velocity 
is compared to the magnitude of the vertical velocity. If the ratio of the terminal 
velocity's magnitude to the vertical velocity’s magnitude exceeds the square root 
of an order of magnitude (i.e., 3.16), the grid point is considered to satisfy 
stratiform conditions. Only regions of the cloud volume where the magnitudes 
of the terminal velocity and vertical velocity exceed 2 and 1 m/s, respectively, are 
examined in order to avoid areas along the cloud boundaries where the 
magnitudes involved are quite small but the ratios large. Above the melting 
layer, only positive values of W are used in the ratio. Finally, below the melting 
layer, in order to capture new convective updrafts that may not yet have 
substantial hydrometeor fallspeeds, any point that has a vertical velocity of at 
least 1 m/s and 0.1 g/kg of cloud water present is considered convective. 

In this study, two criteria that have been adopted in the GCE method are 
applied at the end of all of the other convective-stratiform partitioning methods 
as well. First, surface rainfall rates are thresholded. The threshold is allowed to 
vary between 10 and 25 mm/h depending on the simulation time. Within the 
first hour, the threshold is allowed to vary between 10 and 15 mm/h, and after 
three hours it is allowed to vary between 20 and 25 mm/h depending on what 
the maximum surface rainrate is at that particular time. Any grid column with a 
surface rainfall rate that exceeds the threshold is automatically made convective 
if it already is not. And secondly, an attempt is made to make the convective and 
stratiform regions coherent. This is accomplished by first making all cloudy 
areas ahead of the gust front convective. Next, any cloudy areas that are in front 
of the last stratiform point (in the direction of storm motion) but behind the back 
edge of the convective region are made stratiform. Then, any points classified as 
stratiform that are in front of the leading edge of the convective region are made 



convective. And finally, if the convective region is 5 grid points or less in width, 
any stratiform points that fall in between are made convective. 

Table 2 lists the characteristics of the various convective-stratiform 
partitioning methods. 

3. Results and Discussion 

Models simulations were made for 720 minutes for both the PRESTORM and 
TOGA COARE cases using both Positive Definite Advection (PDA) and fourth- 
order advection (4th) in each case for a total of four experiments. Figure 1 shows 
the vertical cross-sections of estimated radar reflectivity at minute 720 for all the 
runs. All the simulations show a classical squall-line structure (e.g., Gamache 
and Houze 1982) with a leading convective edge, trailing stratiform region with a 
bright band, and new cell growth just ahead of the leading edge as the system 
propagates forward. Reflectivities exceed 50 dBZ and regularly reach 60 dBZ in 
the convective cores of the PRESTORM case for both types of advection. in terms 
of the stratiform region, the 4th order run actually produces on average a slightly 
broader area of stratiform rain. However, the most notable difference between 
the 4th order and PDA can be seen in the variability of the reflectivity field in the 
stratiform region. Not only does the PDA run appear smoother, but the 
reflectivity features are also more coherent. The 4th order run shows a lot of 
small scale variation not present in the PDA run 1 . 

Looking at the TOGA COARE case, many of these same characteristics are 
evident, though there are some differences. Core reflectivities do frequently 
exceed 50 dBZ but never 60 dBZ, and the 50 dBZ values are for the most part 
restricted to below 3 km, unlike in the midlatitude case where they can reach 
upwards of 6 km. Again, the 4th order run produces a broader area of stratiform 
rainfall, but there are also significant differences in the convective region. The 
4th order run tends to produce numerous shallow convective cores that are 


1 As will be seen later in the TOGA COARE hydrometeor profiles, 4th order 
advection can lead to vertical oscillations in the hydrometeor fields, namely 
graupel and snow, not present with PDA. This was also found in another 
current GCE modeling study by Johnson et al. (2001). 



narrow and spaced close together, while the PDA run generates fewer cells that 
are both broader and more likely to penetrate further above the freezing level. 
The sensitivity to advection in the TOGA COARE case not seen in PRESTORM 
may be due to the weaker instability in the tropics that would allow subtle 
differences in advection to have more of an effect. The PDA results agree more 
favorably with the observed Doppler derived reflectivity and vertical velocity 
structures (Jorgensen et al. 1997). Compared to the PRESTORM case, the 
transition from convective to stratiform rain does not appear to be as definitive. 
This is due to the succession of numerous small cells in the tropical case that do 
not significantly penetrate the freezing level, especially with 4th order advection. 
In contrast, the leading convective cells in the PRESTORM case penetrate well 
above the freezing where hydrometeors are then more likely to be carried 
rearward to form a more discerable bright band. 

Time-domain cross-sections of surface rainfall, or Hov-Muller diagrams, 
are presented in Figure 2 for all four simulations. All runs show a change in 
forward propagation speed as the systems mature. The system almost discretely 
jumps to a faster propagation speed at 180 minutes for the PRESTORM case with 
4th order advection. The evolution to a faster propagation speed is similar for 
the PRESTORM PDA run, but the transition is smoother. The shift in 
propagation speed is also accompanied by a broadening of the stratiform region 
(more evident in the 4th order run). The systems are likely transitioning from 
an erect updraft profile to an upshear profile (Ferrier et al. 1996). As for the 
TOGA COARE runs, there is a slight forward acceleration evident after about 150 
minutes in the 4th order run while the PDA run shows a much sharper increase 
at 270 minutes. Also, the system propagates approximately 0.58 m/s faster using 
PDA. The forward propagation speeds were nearly identical for both 4th order 
and PDA for the PRESTORM case. 

3.1 Rainfall Statistics 

Table 3 shows the total amount of rainfall per grid point accumulated over 
all 12 hours of simulation time for all four runs. The totals are similar for each 
case. In the PRESTORM case, PDA produces 1.9% less rainfall than 4th order 
advection, while PDA produces just 1.4% less rainfall in the TOGA COARE case. 
Figure 3a shows time series of instantaneous grid averaged total and stratiform 



rainfall for both PDA and 4th order runs for the PRESTORM case. The stratiform 
amount is obtained via the GCE separation method. The stratiform amounts 
match closely and slowly increase throughout the simulation. Total amounts 
are also quite similar, except for a few distinct episodes of increased rainfall near 
200 and just before 360 minutes in the 4th order run and just after 480 minutes in 
the PDA simulation. Looking at the Hov-Muller diagrams in Fig. 2, these 
episodes can usually be traced to periods when cells develop ahead of the main 
line and then merge with the system. This merging can lead to enhanced 
rainfall (Byers and Braham 1949; Simpson 1980; Wescott 1984; Tao and Simpson 
1984, 1989). 

Figure 3b shows the same instantaneous grid point averaged total and 
stratiform rainfall traces but for the TOGA COARE case. The 4th order run has a 
period of higher rainfall from 120 to 240 minutes followed by a period of higher 
rainfall in the PDA run. The PDA run also has more rainfall in the stratiform 
region from about minute 300 to minute 420 which can be attributed to a few 
cells that persist with greater intensity behind the leading edge in the PDA run. 
Stratiform rain increases more rapidly than in the PRESTORM case then levels 
off. 


Results from the various convective-stratiform separation techniques are 
also shown in Table 3. For the PRESTORM case, the C&H method has the most 
stratiform rain with the GGE and CA20 methods close behind. The Xu method 
then detects a modest amount, while the Vt-W and Steiner methods find very 
small amounts of stratiform rain. The results do not vary significantly relative 
to the type of advective scheme used. The GCE, C&H, and CA20 methods all 
give pretty similar results for both types of advection. These range from 16.9 to 
21.5% stratiform rain inclusive of both advection types. Johnson and Hamilton 
(1988) reported an average stratiform rain amount of 29% for this PRESTORM 
case using surface rain gauge data. As the system matures, the stratiform 
percentage increases both for the model results and the observations. If the 
stratiform rain percentage is computed over the final 240 minutes of the 
simulation, during the mature stage, for the PDA run, the GCE, C&H, Xu, and 
CA20 estimated percentages increase around 10% while Steiner and Vt-W 
increase only 2-5%. Johnson and Hamilton (1988) estimated the observed 
mature stage stratiform component to be between 30 and 40%. This makes all of 



the methods appear too low in their estimates of stratiform percentage though 
the GCE, C&H and CA20 values are reasonably close especially for the mature 
stage. However, Johnson and Hamilton (1988) simply assigned rainrates in 
excess of 6 mm/h as convective in their analysis of the mesonet rain gauge data. 
Applying this same criteria to the model simulations results in a stratiform rain 
percentage of just 4.2% for the PDA run, which is now in very close agreement 
with the Steiner and Vt-W estimates. Obviously the rainfall distribution in the 
model does not match the observations. Sui et al. (1998) found similar problems 
with the model versus observed rainfall. This makes evaluating the relative 
performance of the separation methods using the model versus observed 
stratiform percentages problematic. 

In the TOGA COARE case, all of the methods are consistent in finding 
higher stratiform amounts in the PDA run (3-8% more). This time it is the GCE, 
C&H, and Xu methods that yield similar amounts of stratiform rain, about 35 to 
42% for the 4th order run and about 42 to 49% for PDA. The CA20 method 
shows around 30% and now ranks lower (less than Xu) in a relative sense 
compared with the PRESTORM case. The Vt-W and Steiner methods still give 
the least amount of stratiform rain, less than 20% for the TOGA COARE 
simulations. Short et al. (1997) reported that 40% of the rainfall that occurred 
during active periods in TOGA COARE was stratiform. Active periods are 
associated more with organized types of convection. They used a texture 
algorithm similar to the Steiner method to separate rain maps derived from 
shipboard radar data into convective and stratiform portions. This makes the 
GCE, C&H, and Xu methods appear favorably close to observations in the TOGA 
COARE case. Of course, the same issue of how well the model replicates the 
observed rainfall distribution remains. 

For both the PRESTORM and TOGA COARE cases, the Steiner method 
consistently produces the least amount of stratiform rain and could appear to 
significantly underestimate the amount of stratiform rain. One factor pertains to 
the grid resolutions in the model simulations. The Steiner method was 
designed for a grid resolution of 2 km, whereas the model has a 1 km resolution 
in the PRESTORM case and a 0.75 km resolution for TOGA COARE. If the 
method is applied to a grid with a higher resolution than what it was originally 
designed for, it will detect small scale deviations and identify them as convective 



thus overestimating the amount of convective rain (Steiner et al. 1995). Steiner 
et al. (1995) recommend applying the technique at a level below the bright band 
but not too low so as to allow for greater radar coverage. They also recommend 
using a lower level in mid-latitudes to account for possible evaporative effects. 
In the model, there is no radar restriction, but nonetheless a level of 1.8 km was 
used in both the PRESTORM and TOGA COARE cases. This relatively low level 
means that evaporation is likely to be small. Nevertheless in its current 
application, in the PRESTORM case, regions that the method identifies as 
convective are always close to areas of at least 40 dBZ (Fig. 1). And, in the TOGA 
COARE case, the difference between the Steiner method and most of the others is 
in the treatment of the cells that are behind the leading edge. These are more 
likely to be identified as convective by the Steiner method and Vt-W methods. 
Occasionally the GCE and Xu methods identify parts of these as convective too 
(Fig. 1). So, based on the simulated reflectivity structures and the possible bias 
towards high rainrates in the model, the Steiner method appears to perform 
reasonably. 

3.2 Area Cloud Coverage 

Vertical profiles of cloud coverage are shown in Fig. 4 for both types of 
advection for each case. Grid points that had total hydrometeor contents in 
excess of 0.01 g/m3 were considered cloudy. In the PRESTORM case, below 7 km 
the profiles are very similar, while at upper levels, on average, the PDA run 
shows a broader area that is associated with the non-surface-precipitating portion 
of the anvil. The difference is largest at 12 km, on the top portion of the anvil. 
The mean width of the cloudy area is 280 km for 4th order advection while PDA 
averages 304 km. In the TOGA COARE case, the situation is similar with the two 
profiles nearly identical below about 5.5 km and PDA broader above that level. 
The average cloud areas in the TOGA COARE case are 142 km for 4th order and 
161 km for PDA. In addition, average cloudtops in the PDA run are a kilometer 
higher. This is likely a result of the taller more vigorous nature of the cells in 
the PDA simulation that allows more condensate to penetrate higher above the 
freezing level. 

Time series of convective and stratiform rain area (not shown) were 
computed for both cases and both advection types using each of the separation 



techniques. In the PRESTORM case, there is an overall trend for the stratiform 
area to slowly but steadily increase in size. The increase flattens out over the 
final 2 hours of simulation time. The average mean width of the stratiform 
region over the 12 hour simulation time for all 6 methods was 47 km for 4th 
order and 43 km for PDA (Table 4). In the PDA run the mean 12 hour stratiform 
width ranged from 50 km for the C&H method to 29 km for Vt-W. Convective 
area increased much more rapidly. The growth tapered off a little more quickly 
in the 4th order run. The resulting mean width of the convective region for the 
12 hour simulation time was 81 km for 4th order and 95 km for PDA. The 
average convective width varied from 87 km for C&H to 110 km for Vt-W in the 
PDA run. Because the coherence approach was applied, the CA20 method 
produced an average convective width of 91 km for the PDA run due to the 
forward anvil. This agrees quite well with the ensemble average. 

The TOGA COARE runs reveal a similar trend in the stratiform area to 
the PRESTORM runs, a slow steady increase over the course of the simulation. 
The resulting average mean 12 hour stratiform width for the 4th order and PDA 
runs (Table 4) was 62 and 56 km, respectively. In the PDA run, the mean 12 hour 
width varied from 43 km for Vt-W up to 67 km for C&H. The convective width 
in the TOGA COARE runs increased more slowly than in the PRESTORM case 
and leveled off much sooner at a much smaller characteristic width (partially 
due to the forward anvil in the PRESTORM case). The 12 hour mean average 
values were 20 km for 4th order and 27 km for PDA for TOGA COARE. So in 
both cases the convective ratio was slightly higher for PDA. The characteristic 20 
km convective width assumed in the CA20 approach yields 12 hour means that 
in the tropical case again agree quite well with the ensemble, 21 and 28 km 
respectively for 4th order and PDA. 

3.3 Cloud Mass Fluxes 

Accumulated mass fluxes for the PRESTORM case are shown in Fig. 5a for 
both types of advection. The fluxes are subdivided into convective and anvil 
regions using the GCE method. The convective profiles are fairly similar 
showing positive at all levels with a single maximum in the 4 to 5 km range. 
Though the maximum is larger for the PDA run, at upper levels, the 4th order 
convective profile is larger. In the anvil region, both schemes show single peaks 



in upward fluxes at 9 km and single downward peaks between 2 and 3 km. Aloft, 
the upward peak is larger for PDA while at low levels the two schemes produce 
similar downward peaks but at slightly different levels. If the convective and 
anvil fluxes are combined however, the total profiles are nearly identical at all 
levels for both types of advection suggesting the differences in the component 
profiles may be attributed to the convective-stratiform partitioning. Time-height 
cross-sections of maximum downdrafts (not shown) reveal that peak updrafts are 
found initially around 9 km during the early stages of the simulations then 
transition down to 4 km over the later stages for both types of advection. 

The mass flux profiles for the two TOGA COARE simulations are shown 
in Fig. 5b. The convective and anvil profiles, again using the GCE separation 
method, are generally similar to those found in PRESTORM. The convective 
peaks are sharper and lower around 2 to 3 km, about the same level as the 
downward peaks in the anvil profiles. If the components are combined 
however, the total accumulated profiles do show some differences. Between 2 
and 6 km, there are more net upward fluxes in the 4th order run. At upper 
levels in the anvil region, the PDA simulation contains significantly more 
upward fluxes. These additional fluxes are associated with cells behind the 
leading edge that are much more vigorous in the PDA run. Time-height cross- 
sections of maximum updrafts (not shown) indicate that after 240 minutes, peak 
updrafts occasionally exceed 3 m/ s in the 4th order run above the freezing level 
but are usually on the order of 1 to 2 m/s. On the contrary, the PDA simulation 
contains updrafts regularly on the order of 3 m/s but sometimes as high as 5 to 7 
m/s above the freezing level. Furthermore, in the PDA run, a few 2 to 3 m/s 
updrafts are able to penetrate above 13 and 14 km, whereas after 240 minutes, 2 
m/s updrafts in the 4th order run are confined below 13 km. 

Figure 6a shows profiles of accumulated mass flux separated into 
convective and anvil components for the PDA PRESTORM run using each of 
the six different separation techniques. All of the profiles are qualitatively 
similar, and the results appear independent of the advection type. The 
convective profiles all show net upward mass flux throughout the depth of the 
troposphere with a single maximum. The anvil profiles all show net upward 
mass flux in the upper troposphere and net downward mass flux in the lower 
troposphere. The heights and magnitudes of the maxima vary significantly from 



technique to technique. The height of the convective maximum varies between 
3.5 and 6 km while the magnitude varies between 110 and 130 g/cm2/grid. The 
stratiform maxima have comparable differences in height but even greater 
variation in magnitude. In general, the techniques behave as expected especially 
at middle and upper levels where the Vt-W method has the convective profile 
with the most upward mass flux and the C&H convective profile has the least. 
The other methods fall in between in the same relative sequence as their 
stratiform areas except for the CA20 method which differs significantly from the 
rest especially at low levels. Based on cross-sections of vertical velocity, there is 
evidence for two main downdrafts behind the leading edge updraft. The first 
downdraft is located about 5-10 km to the rear of the leading edge updraft, and 
the second is generally about 20-25 km behind. All the methods include the first 
downdraft nearest the leading edge updraft as part of the convective region. The 
convective area of the CA20 method does not reach to the second downdraft 
while the other methods occasionally include it as part of their convective 
regime. This accounts for the larger negative mass fluxes in the low level anvil 
profiles of the CA20 method compared to the others. With the second main area 
of downdrafts effectively in the anvil region, the low level convective region of 
the CA20 method shows more net upward fluxes resulting in a maximum at 
lower levels. Some of the methods include some very strong updrafts aloft as 
part of the upward mass flux in the their anvil region [see Section 3.6]. 

Compared to the PRESTORM profiles, the TOGA COARE convective and 
anvil mass flux profiles generated by the various separation techniques (Fig. 6b) 
are less variable at low levels but more variable aloft. Again, the results are 
consistent for both types of advection. Almost all of the methods place the 
convective peak upward mass flux at a level of 2.5 km. The peak convective 
magnitudes vary from 40 to 65 g/cm2/grid. Likewise the low-level downward 
fluxes in the anvil profiles are also quite similar. The GCE, C&H, and CA20 
methods have nearly identical convective and anvil profiles throughout the 
depth of the system. The real variation between the methods appears aloft, 
above the freezing level. The Xu, Steiner, and Vt-W methods increasingly 
categorize upward mass fluxes above the melting level as convective whereas 
the other three have slightly negative net convective mass flux above the 
melting level. The Xu method is based on vertical velocities. Occasionally 
vertical velocities near the freezing level and to the rear of the leading edge were 
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sufficient to satisfy the convective requirements stipulated in the Xu method. 
These updrafts were associated with cells behind the leading edge. While the Xu 
method categorized parts of these cells as convective from time to time, the 
Steiner method was moie likely to identify them as convective based on their 
reflectivity structure. Frequently these cells contained updrafts aloft that were on 
the order of 1 to 2 m/s. Consequently, the Vt-W method was most likely to 
categorize these cells as convective leading to a fairly large diversity among the 
profiles above the freezing level including changes in shape. 

3.4 Apparent Heating (Qj) and Moistening (Q 2 ) 

Q1 profiles between 4th order and PDA in the PRESTORM case are nearly 
identical. Figure 7a shows that the total heating matches over the complete 
depth of the troposphere for the two different advection types. The total heating 
shows a double maxima with a mid-level maximum at 7 km and a smaller 
secondary one near 13 km. There is heating throughout the depth of the 
troposphere but cooling above the tropopause near 15 km. The convective and 
anvil heating profiles, partitioned according to the GCE method, are also almost 
identical. The convective profiles show predominantly heating throughout with 
a main maximum near 5.5 km. A relatively weak secondary maximum is 
evident near 12 km, and the cooling above the tropopause is minimal. The 
anvil profiles show cooling below 5 km, though the freezing level is near 4 km, 
and heating aloft with the largest maximum near 9 km and a secondary one 
nearly as large near 14 km. The cooling above the tropopause is more 
pronounced than in the convective profiles. All of the profiles have been 
normalized with respect to total surface rainfall. 

Convective and anvil profiles for PDA are broken down into their 
individual components in Fig. 8a&b. The convective components are quite 
similar for the two types of advection. Condensation is the dominant term with 
a maximum near 4 km, close to the freezing level. Above this level, 
condensation drops off while deposition, the second largest term, increases to a 
maximum near 9 km. Sublimation is relatively small and peaks at the same 
level as deposition. Evaporation increases towards the ground and becomes 
larger than condensation only near the surface. Freezing and melting are 
combined into one curve and are of equal magnitude on either side of the 



freezing level. The weak secondary maximum in heating near 12 km is 
attributed to eddy flux convergence. This term is also responsible for the cooling 
above the tropopause. It is also a relatively important term near the freezing 
level and again near the surface. Radiation effects are negligible in the 
convective region. 

The anvil components, shown in Fig. 8b, are also similar for both types of 
advection except around 9 km. Condensation is still significant, but deposition 
and sublimation are now the dominant terms. Condensation is fairly uniform 
from about 1 to 8 km. Evaporation is quite weak above 4 km but is the dominant 
term below 3.5 km. Deposition dominates the heating above 6 km, but it is 
nearly compensated for by sublimation. The peak values in deposition and 
sublimation are much greater in the PDA run. However, the net difference 
between the two is the same as the 4th order run resulting in the same net 
heating. Freezing is negligible while melting is significant over a 2 km depth 
below the melting level. The eddy flux term is significant below the melting 
level, but is quite large around the tropopause. It alone is responsible for the 
elevated secondary peak in anvil heating and also for the large cooling seen 
above the tropopause. The radiation term does show cloud top cooling and cloud 
base warming in association with the anvil, but it is still in comparison a 
relatively minor term in the anvil region. 

Unlike the PRESTORM case, the TOGA COARE Ql profiles for the two 
advection schemes do show some differences (Fig. 7b). Again using the GCE 
separation method, the convective profiles have similar shapes but the 4th order 
profile has a larger magnitude of peak heating around 3 km. The anvil profiles 
are similar below 10 km, but above that level the PDA profile shows significantly 
more heating. These differences are reflected in the total heating profiles. 
Compared to the PRESTORM profiles, the TOGA COARE convective profiles 
show a similar magnitude in peak heating when normalized with respect to 
surface rainfall, but the elevation of the peak is much lower near 3.0 km 
compared to 5.5 km in the PRESTORM case. Also, whereas the PRESTORM 
convective profiles were nearly all heating, the TOGA COARE profiles do show 
some significant cooling from 7 to 12 km above the lower level heating. But 
unlike the PRESTORM case, there is no strong heating or cooling associated with 
the tropopause. The anvil profiles show strong cooling, stronger than in the 



PRESTORM case, from below 4 km down to almost the surface. The peak 
heating aloft is also stronger than in PRESTORM though at about the same level, 

7-9 km. This heating instead shows up in the convective profiles for 
PRESTORM. 

The convective components for the PDA TOGA COARE profiles, shown 
in Fig. 8c, are similar to the PRESTORM components in that condensation is the 
dominant term. However, in the tropical case this term is so dominant that 
most of the other terms are practically insignificant by comparison. With a 
higher freezing level and weaker instability in the tropics, deposition and hence 
sublimation become far less important. The same is true for freezing and 
melting. The moister environment also helps to reduce evaporation in the 
convective region though it is non-negligible. The eddy flux term is smaller too 
compared to the midlatitude system. As with the PRESTORM case, the 
convective components are independent of the advection. 


Figure 8d shows the anvil Ql components for the PDA TOGA COARE 
simulation. Heating is almost evenly divided between condensation below and 
deposition aloft while sublimation is far smaller than in the PRESTORM case. 
Eddy fluxes are slightly weaker than for PRESTORM, but evaporation is quite 
large below 5 km. Above 9 km, deposition and to a lesser degree sublimation are 
greater with PDA. Radiational cloud-base warming is also larger and shifted 
upward by 1 km in the PDA simulation. Likewise eddy flux warming just above 
cloud top is stronger and 1 km higher. This is consistent with the taller cells in 
the PDA simulation for TOGA COARE. 

Normalized convective and stratiform Ql profiles from the various 
separation techniques are shown in Fig. 9a&b for the PRESTORM case. Again, 
the results are similar for both types of advection. The overall shapes of the 
profiles are very similar between the methods with the only real differences 
being in the magnitudes at mid-levels. Mid-levels are where the greatest 
variation exists; Steiner and Vt-W have the largest convective heating (and 
smallest anvil heating) while C&H has the least. The Xu, GCE, and CA20 
methods fall in between in descending order of convective heating. The biggest 
differences are at 7 km in the 4th order run where convective heating ranges 
from 0.75 (deg /mm) for the Steiner and Vt-W methods down to 0.43 for C&H. 



Below the melting level and above 9 km, the methods are all very close with 
maximum differences on the order of 0.1 (deg/mm). 
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Results from the two TOGA COARE simulations are shown in Fig. 9c&d. 
The separation techniques show very close agreement below the melting level. 
The narrowness of the results in this region can be attributed to the lack of any 
strong microphysically induced heating or cooling beyond the vicinity of the 
leading convective edge. Since all of the methods assign the region in and 
around the leading edge as convective, there is very little variation in the results 
even though the convective areas may differ significantly. Unlike the 
PRESTORM case where strong pockets of cooling, due to evaporation or even 
melting, can exist well behind the leading edge below the melting level, the 
moist environment of the tropics mitigates this effect. Hence the Steiner and Vt- 
W methods, which on average have significantly larger convective regions, yield 
only slightly less cooling in the anvil region compared to the other methods 
below the melting level. The results are different at middle and upper levels. 
The GCE, CA20 and C&H methods all produce similar convective and anvil 
heating above the freezing level. The Xu method differs slightly from these 
assigning more of the heating aloft to the convective region. The Steiner 
method is next with even more of the heating as convective, while Vt-W has the 
most heating aloft in the convective profiles noticeably more than even Steiner 
This is different from the PRESTORM case where the two methods are nearly 
identical at all levels. Unlike below the melting level, there are regions of 
significant heating aloft to the rear of the leading convective edge where there 

are pockets of condensational and depositional heating associated with old cells 
to the rear of the leading edge. 

Comparing the Q2 profiles for the two PRESTORM runs (Fig. 10a) shows 
that total moistening for 4th order and PDA is nearly identical. The convective 
and anvil profiles are also very similar. The TOGA COARE Q2 profiles (Fig. 10b) 
are also quite similar between 4th order and PDA. The anvil profiles are a very 
close match while 4th order has slightly more drying between 2 and 4 km in the 
convective profile and hence the total profile. Both of the TOGA COARE anvil 
profiles show moistening below 4.5 km, similar to the PRESTORM profiles only 
stronger. Both the TOGA COARE and PRESTORM convective profiles show 
drying mainly below 4.5 km with similar peak magnitudes. The TOGA COARE 



profiles peak near 3.5 km while the PRESTORM profiles peak closer to 2 km. 


The convective and anvil Q2 profiles show greater variation in the 
PRESTORM case than in the TOGA COARE case (Fig. 11). This is again due to 
the lack of strong microphysics behind the leading edge below the melting level 
in the TOGA COARE case where Q2 has its largest amplitude. Most of the 
methods are in close agreement in the PRESTORM case except for CA20. The 
CA20 method shows the largest deviation from the other methods with the 
strongest anvil moistening and convective drying for both runs. All six 
separation techniques are similar to each other in the two TOGA COARE 
simulations below 4 km differ slightly at middle levels, however. 

3.5 Hydrometeor Profiles 

Total hydrometeor profiles for both forms of advection and both cases are 
shown in Fig. 12. The two PRESTORM runs show very close agreement for all 
five hydrometeor species. The snow profile for 4th order advection has a slightly 
larger peak value. The TOGA COARE runs show more variation. At low levels, 
the 4th order rain and cloud profiles are slightly larger though the shapes are 
similar. Mid-level graupel contents are higher in the 4th order run, but aloft the 
PDA snow, graupel and cloud ice contents are substantially larger which is 
consistent with the deeper cell structure. Also, the PDA graupel and snow 
profiles are much smoother. 

The GCE method is the only one that directly utilizes hydrometeor data to 
help differentiate between convective and stratiform areas. Though difficult to 
obtain in real world applications, it can be useful for identifying convective 
regions. Large values of cloud water or cloud ice are associated with significant 
updrafts that are identified with convection. Total hydrometeor profiles for the 
convective and anvil regions were computed for each technique (not shown). 
Differences in cloud water content between the various methods for the 
PRESTORM case are not very large. The Steiner and Vt-W methods have the 
most convective cloud water (97%) while C&H has the least (74%) The GCE 
method is in the middle (84%). In all of the methods, the percentage of cloud 
water that is convective is high. Aloft, cloud ice is roughly equally divided 
between convective and anvil. Values range from 83% convective for Vt-W to 



only 47 /o convective for C&H. Another indication of convection is the presence 
of hail. For the Vt-W and Steiner methods, around 85% of the hail falls within 
the convective region; for Xu it is 68%, for the GCE and CA20 methods, about 

58 /o, and for C&H, only 50%. The behavior between PDA and 4th order 
advection is the same. 


In the case of TOGA COARE, nearly all of the cloud water below 3 km (2 
km for PDA) is found to be convective in all the methods. Cloud water is 
present at mid-levels in both the PRESTORM and the TOGA COARE case, but 
unlike the PRESTORM case, there is a lot of variation between the techniques 
regarding the mid-level cloud water. This is due to the the presence of 
significant amounts of cloud water further behind the leading edge in TOGA 
COARE. Cloud water ranges from just 48% convective in the C&H method up to 
88% for Steiner. GCE has around 60%. The GCE, C&H, and CA20 methods have 
essentially no convective cloud ice (less than 7%) with Xu only a small amount, 
about 12%. The Steiner method has a fair amount with 37% while the Vt-W 

method has the most at 47%. The values are for the 4th order run; the PDA 
results are similar. 


3.6 


Contoured Frequency with Altitude Diagrams (CFADs) 


CFADs are frequency distributions at each altitude summed together as a 
function of height to make a statistical diagram of a particular field (Yuter and 
Houze 1995b). Steiner et al. (1995) constructed CFADs of radar reflectivity from 
Darwin, Australia, for convective, anvil, and total storm volumes for the entire 
month of February 1998 to verify their convective-stratiform separation scheme 
(i.e., the Sterner method). A narrower distribution with peaks at lower 
magnitudes (especially at low levels) is expected for the stratiform region, which 
they found. They also tested their separation algorithm by computing CFADs on 
data from a CaPE (Convection and Precipitation/Electrification Experiment) case 
where doppler derived vertical velocities were available. They found that the 
velocities were consistent in showing a much narrower CFAD for the stratiform 
region and a broad convective CFAD having a significant portion of the vertical 
velocities greater than +-5 m/s. The corresponding reflectivity CFADs were also 
compatible with those generated for the Darwin data. 



CFADs were constructed for both the PRESTORM and TOGA COARE cases 
for each of the various separation algorithms. The PRESTORM results indicate 
that each of the algorithms was able to produce convective and stratiform 
reflectivity CFADs that show larger peak modes in the convective CFADs, which 
is consistent with the results Steiner et al. (1995) used to verify their algorithm. 
CFADs for the Steiner and GCE methods are shown in Fig. 13a-d for PDA. At the 
lowest levels, the peak convective modes were close to 50 dBZ for GCE, C&H, 
CA20 and Xu down to 45 dBZ for Steiner and Vt-W for PDA. The stratiform 
modes range from 20 to 25 dBZ for GCE, C&H, CA20 and Xu to 20 dBZ for Steiner 
and 15 to 20 dBZ for Vt-W also for PDA. In addition, all of the stratiform CFADs 
from the PRESTORM case showed evidence both of a bright band and of 
decreasing magnitudes below the bright band indicative of evaporation. The 
effects of low-level evaporation are not seen in any of the convective CFADs. 
One distinct difference from the results of Steiner et al. (1995) is that most of the 
stratiform CFADs generated by the C-S techniques in this study are broader than 

those of the convective region. The Steiner and Vt-W methods are the 
exception. 

Analysis of the reflectivity CFADs for the TOGA COARE case reveals 
similar trends to those from PRESTORM. Results from Steiner and GCE are 
shown in Fig. 13e-h for PDA. The GCE, C&H, and Xu methods are in close 
agreement all having convective peak modes just under 50 dBZ for both types of 
advection and peak stratiform modes near 20 dBZ for 4th order and near 30 dBZ 
for PDA. The CA20 method, however, is now a bit different. Though it still 
agrees with the others in terms of stratiform distribution peaks, it now has peak 
convective modes closer to 40 dBZ. Furthermore, the convective CFADs for the 
CA20 method are now much broader than those of GCE, C&H, and Xu. These 
results are likely attributable to the difference in cell structure between the 
PRESTORM and TOGA COARE cases. There were more numerous cells behind 
the leading edge in the TOGA COARE simulations. The Steiner method has 
smaller convective and stratiform peak modes than the first four methods, about 
37 for the convective peaks and 18 and 25 dBZ for the 4th order and PDA 
stratiform peaks respectively. Finally, the Vt-W method shows the smallest 
convective peaks near 35 dBZ and fairly small stratiform peaks from 20 to 22 
dBZ. Thus all the methods are again able to produce larger peak values for the 
convective region. All of the stratiform CFADs do show evidence of low-level 



evaporation for the 4th order simulation, yet none of them show it for PDA. 
Only the CA20, Steiner and Vt-W methods show convective distributions that 
are close to being as broad as the stratiform distributions. 


Vertical velocity CFADs for the PRESTORM case all have broader 
distributions in the convective region though the distributions are overall much 
narrower than those in Steiner et al. (1995). The vertical velocity CFADs for Vt- 
W and GCE for PDA are shown in Fig. 14a-d. At +5 m/s, the contour frequency is 
close to 1 /o at most levels in the current study whereas Steiner et al. (1995) has 
the frequency exceeding 10% at mid-levels. The reason for the discrepancy is that 
in the current study, clear areas, especially ahead of the leading edge above new 
convection, are counted as part of the convective region and skew the 
distributions closer to zero. Also, the horizontal domain in Steiner et al. (1995) is 
much smaller. The vertical velocity CFADs reveal that some of the separation 
techniques contain a few occurrences of strong vertical velocities in the 
stratiform region. These include a few 15 to 20 m/s updrafts in GCE, Xu, and 
CA20, and some that reach 25 to 30 m/s in C&H. The Steiner method has some 
10 to 15 m/s updrafts in the stratiform region, while the Vt-YV method contains 
no updrafts over 5 m/s in the stratiform region. These outliers are infrequent 
and do not affect the overall results, though they are clearly misclassified. Xu 
(1995), however, states that significant vertical velocities above the melting layer 
not associated with convective rain rates should be classified as stratiform as they 
are predominantly in the upper troposphere. The largest outliers do in fact occur 
about the freezing level, but they do not satisfy the condition Iwl « Vice put 
forth by Houze (1993) and Steiner et al. (1995) where Vice is of the order 1 to 3 
m/s. The broadest convective and narrowest stratiform CFADs are those from 
Vt-W. Steiner is next, followed by the others which are similar to each other. 
The results are the same for both types of advection. 

CFADs of vertical velocities for the TOGA COARE case, though slightly 
narrower, show similar characteristics compared to PRESTORM. Results for the 
GCE and Vt-W methods are shown in Fig. 14e-h for PDA. Convective CFADs for 
the 4th order run (not shown) do show an enhanced region of positive vertical 
velocities near 2 km compared to those for PDA. This may be a result of the 
more numerous shallow cells. Again, the Steiner and Vt-W methods have the 
broadest convective CFADs and the narrowest stratiform ones. 
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3.7 Microphysical Processes 

In addition to the true CFADs of vertical velocity and radar reflectivity 
another type of contoured diagram was made using microphysical data from the 
model. These diagrams, termed CDADs for Cumulative Distribution with 
Altitude Diagrams, were constructed by partitioning the microphysical processes 
of condensation, evaporation, deposition, sublimation, freezing and meltinv 
according to vertical velocity bins (the same used in the vertical velocity CFADs° 
again as a function of height, cumulatively. Finally, each bin was normalized bv 
t e total cumulative value for that process independent of the region (i e not 
convective or stratiform but total) involved so that the diagrams show as a 
unction of W and Z the distribution of where a particular process occurs. These 
lagrams were constructed for the entire domain and for the convective and 
stratiform regions for each separation technique. Results for the entire domain 
are shown in Fig. 15. The total accumulated distributions of each the 
nt.crophys.cal processes are also partitioned into convective and stratiform 
fractions for each separation technique. The results are listed in Tables 5 6 7 


In the PRESTORM simulations, the bulk of the condensation takes place 
in updrafts between 2 and up to 14 m/s between 1.5 and 6 km in altitude. All of 
the convective-stratiform separation techniques classify the bulk of the 
condensation as convective ranging from just over 70% for C&H in the 4th order 
run to over 97% for Steiner and Vt-W (see Tables 5&6). Based on the convective 
and stratiform CDADs (not shown), all of the methods show that evaporation 
can e separated into two rather distinct regimes. The first regime is the 
convective regime wherein evaporation occurs from the surface to up over 6 km 
with a significant portion from downdrafts stronger than -3 m/s. The convective 
distribution is skewed to the downdraft side extending to regions over -5 m/s 
This evaporation is attributed to downdrafts along the periphery of rising updraft 
cores. The stratiform regime, on the otherhand, is narrower and confined below 
km. The distribution is centered on the zero line and becomes even narrower 
near the surface. Evaporation in the stratiform region is dominated by rain 

“ "If* ’ en ”‘ h the anvil ' ° vera11 for the PRESTORM case, CCE, C&H, Xu and 
0 divide evaporation roughly equally between convective and stratiform. 



Steiner and Vt-W place three-quarters of all evaporation in the convective 
region Tables 5&6). The range being from 44% convective for C&H up to 81% 
for Vt-W inclusive of both advection types. 

The results for deposition are more variable in the PRESTORM case The 
total distribution lies mainly between 5 and 12 km and is skewed towards weaker 
up rafts though a significant portion extends past 10 m/s. The portion below 
a out 5 m/s pertains to the stratiform region, while the convective portion of the 
distribution is centered between 5 and 7.5 m/s. This accounts for the shape of the 
overa istn ution. The partitioning of deposition into convective and 

S 13 ^,v, ,n ranses from i ust 0VCI 30% convective for C&H in the PDA run to 
over % convective for Vt-W in the 4th order run. Overall, the methods treat 
eposition as half convective and half stratiform. There is a tendency for 
epos, t, on to be more convective with 4th order advection. The difference is 
» on average. Sublimation predominantly occurs between 4 and 11 km in 
owndra ,s of less than 5 m/s with the distribution skewed towards weaker 
s own ra (v xcept for the Steiner and Vt-W methods, 75-80% of sublimation is 
stratiform. Steiner and Vt-W are slightly lower between 60 and 70%. 

Freezing is treated as a convective process by all of the convective- 

and VtW "h 7 qUeS ; 11 15 almOS ‘ com P ,etel y “nvective (over 96%) in Steiner 
. "7 r .“ nS y convectlve in *e rest with C&H the lowest at 65% for 4th 

, in e istribution is slightly different between 4th order and PDA The 
distribution extends from the freezing level near 4 km up to 9 km The 
, erence between the 4th order (not shown) and PDA distributions comes from 
he contribution by stronger updrafts. For updrafts over 7-8 m/s, most of the 
Zing occurs within a kilometer of the freezing level in PDA whereas it s 
spread more uniformly over 3 to 4 km with 4th order advection. The results 
vary in the assessment of melting. In C&H, melting is nearly half convective 

, '" er a " d V »- W “ i5 ° Ver 80% convective. Convective melting has a 

wider d,s ribution from -10 to 10 m/s while melting in the stratiform region is 

n me etween -5 and 4 m/s. The stratiform distribution peaks near 0 m/s 

* 7 ™ nvectIve distribution is spread from 0 to -5 m/s. Overall, most of the 

Z “ C ° nVeCHVe Ca,e8 ° ry ,h ° U8h 3 occurs in 
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With TOGA COARE, the results for condensation are similar in terms of 
the partitioning but differ in terms of distribution. Again, condensation is found 
to be predominantly convective. Vt-W and Steiner are the most convective (up 
to 92% for Vt-W in the 4th order run). The C&H method considers only 55% of 
the condensation to be convective in the PDA run. In general, the percentage of 
condensation that is found to be convective is at least 10% lower per method 

PRESTOPM “ STORM <=“ The distrib “«°n Offers significantly from the 
PRESTORM results. Condensation has two peaks at lower levels that merge into 

a single peak at mid-levels. This structure is similar for both 4th order and PDA 
The convective distributions show an ascending axis that extends from 10 m/s at 
m up to 1 m/s at 4 km. The stratiform portion of the condensation 
distribution is narrow and lies along the zero axis below 2 km. It widens above 2 
km and extends all the way up to 8 km with a peak just above 4 km at 1 m/s. The 
separation methods produce consistent convective and stratiform patterns of 
condensation as functions of height and updraft. There is considerable variation 
in the partitioning of evaporation between the methods. The GCE, C&H, and Xu 
methods lean heavily towards the stratiform side with percentages over 70% 
stratiform. Steiner and Vt-W are slightly convective with percentages around 
convective. CA20 falls in between with near 40% convective The 
distributions for the convective, stratiform and total regions are similar 
Evaporation extends from the surface up to 6 or 8 km (slightly higher for 
strati orm). It is centered on 0 m/s and falls between 3 m/s with the stratiform 

distribution slightly narrower. Again, the structures are similar among the 
various methods. & 


As in the PRESTORM case, the characterization of deposition as being 
either convective or stratiform varies widely among the separation methods 
The differences for TOGA COARE are even larger with C&H having as little as 
3.7 /„ convective (PDA) and Vt-W up to 67.5% convective (PDA) The 
distribution lies mostly between 5 and 13 km and is skewed towards zero with 
significant values below 8-9 m/s. Among the methods that have a meaningful 
convective amount, the convective side is broader while the stratiform portion 
IS taller and narrower. Sublimation likewise shows a tremendous variation in 

48V°Xrw7 / Th at H 0rm K deSi8nati0n '' aryin8 fr ° m 2 ' 6% COnVective (C&H) to 
( ‘~ W> - The dlstribu tion is centered on zero m/s and varies between 5 

m/s. 



variation"'!!" “ '° n8er P red <™™"«y convective. The 

now considerable between the methods. The most significant 

changes occurred with the C&H and GCE methods. With the GCE method 2 

T^TroATOV° s !!”LtTlHd m 8ai ;v n PRESTORM (4,h) to 24% in 

rr p — ” " 

tot " fe br ° ad t 6 ’ USt abOW the free2i " 8 ,eVel atlribUted 

that the „ 8 . Th convective portion is rather shallow, however so 

toward t 1 "r t ; and ieans 

PRESTORM and TOGA COARF 8 • , b 8 H §e between the 

. . ■ COARE cases involves the melting All of the 

techniques agree in substantially lowering the Derrenu,™ t § 

in the TOGA COARF - j . 8 P ge of convective melting 

(PDA) T1? A ! T ]USt ° Ver 3% f ° r C&H < 4th ) *> *1% ^r Vt-W 

v )• e greatest reduction occurs in the C&-H C'cv: v , ^ 

and is on the order of 40% or more. Stemetand V t ' w“ a" t me ‘ h ° dS 

reduced, near 35% and 20 to 30% respectively This -hift- „• T° n ° >KeaDly 

convective process in PRFSTORM t meltan 8 from °< a 

distribution of melting is more confined in TOGA COARF 7 

PRESTORM with respect to updraft and vertical extent 7eltl! ^7'° 

1.5 km of the melting level. In the 4th order run he d slu ^ "“h 1 " 

o m/s and spans 5 m/s The PDA i ' tb e distribution is centered on 

the convective regin Sh ° WS 3 UP * 10 m/s «»* “mes from 
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Latent Heating Retrieval 


enlLlriTltt^T/o' Id d d° U C d SyS ‘ emS ^ ‘ ht ' ^ 

Heating algorithm (Tao et al. 990- 1 993 . 2000 ) r J ° ddard Convective-Stratiform 
observed convective systems bv u ^ hea,in « P rofii “ for 

observed stratiform" ” 

normalized convective and stratiform heallg 7 f e rZ7T 

advantage of the characteristic shapes in the vertLl ofl7 t T, ^ 
convective and stratiform • he Iatent heating in the 

gions of mesoscale convective systems. Tao et al. 



needed to te accurate to within io% 
ate retrieval. This is based on the fact that the level of maximum 

S u e f fi“n^^ d X W me\ttllT Sin8 Sfratif0rm fraC,1 ° n and 10% is 
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The sensitivity of the retrieved heating profiles to the choirp nf 
technique can also be tested and is shown in Fig. 16 for the PRESTORM ^ 
TOG A COARB PDA cases. Treating the model sibilations Is suTsI^ fort 
real system and using the observed stratiform amount of t 9 "/ n h ? 

Hamilton 1988) for the PRFSTORM r* , f ^ ° ^ ohnson and 

heating profile in the model is although thetcT C^JTcAM "'dTl 
profiles are very good approximations The VhW and S te f ^ 

particular method was in the simulation to the chosen stratiform T ' 

— 1 r ** - ^srz : 

nr t “ ,y - ™ 

c r'ts a " d y r ValUe ° f 40% StraHform < Sh ° rt <* al. 199/1, the GCE 

CA20 aid vl-bblmlsbreTllttlf tooTt^f “ “* ' rUe m ° del h^" 8 ' ^ 
Steiner again produces an exaggerated ^ 

r:bCedbr:;r rr — z 

are going to produ ^ ^ ^ ^ 

methods may yield sienificantlv Jiff , , 8 different separation 

y yie m significantly different stratiform percental ,> ; c o 

normalized profiles that matfpr c , P centa 8 es ' lt )S the ram 

and cooling below and a large stra:“,bmbrb“ at,n8 ^ 

=“s 

- - - - — r zsvj" ■ sssr,™ r r 
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ore so t e teiner stratiform profiles significantly diverge from the other 
ra t orm pro lies. Their shapes appear amplified with respect to the others 
IS IS a direct result of having small stratiform rain amounts by which to 
normalize by in combination with significant heating and cooling present in the 
simulated anvil region with little or no associated surface precipitation. This can 
ave a arge impact on the retrieved profiles. The TOGA COARE profiles also 
s ow strati orm profiles that appear to be amplified relative to the rest This 

I™ 20 “ d Steiner P rofil « are involved. The V.-W profiles, however, 

. 3 . ere can a,so be fundamental changes in the basic shapes not related 

to just division by a smaller denominator. 

4.0 Summary and Conclusions 

* 1 d 1 en * . C ° n V6C * 1 Ve ' S tra t] for m separation techniques were compared and 

uae using 2D numerical simulations of a tropical (TOGA COARE) and 

Zt ::; C °“ (PRESTORM, *>““ Six techniques inlL 

(1995) a y Tr * 0P: ChUrChlU and Houze < 1984 >- Tao et al. (1993b), Xu 

Vt-W method b T 3 "' 7 (e ' 8 ' Ca " iaUX '" 4) ' S ' einer et ah < 1995 >- a ™ a new 
method based on the premise that the terminal velocity of hydrometeors 

el a| Se i995- H ‘° “ ” Stn,tif ° rm P reci P itati » (Houghton 1968; Steiner 

' . ' ° UZe 1997) The simulations were made using two different 

numerical advection schemes: 4th order and PDA. The overall Active was to 

ssess the sensitivity of partitioned quantities such as rainfall, mass flux Q1 Qa 
rn.crophys.es and latent heating retrieval to the use of different convert!' 
separation techniques and to determine the effects of the numerical! 
i. cherries on the simulated systems. 

effect ' erms ,° f rainfal1 ' il was found tbat overall rainfall was not significantlv 
effected (less than 2%) by the numerical advection for either case. However 

esimae stratiform percentages varied significantly among the separation 
cThZ'gcT 4 'b T" " PRESTORM and 12 to “ 9% m TOGA COARE. The 

Steil al V. wTh ? '“It' 0 Pr<>dUCe ‘ he WgheSt Strat,f0m ’ P° rt '° n and 

and Vt-W the lowest. Johnson and Hamilton (1988) reported a stratiform 
rain amount of 29% for the PRESTORM case using gauge data making a„ 7 
techniques appear too low in the model. However, applying the convective 
ramrate threshold used by Johnson and Hamilton (1988) in the model results in 



a stratiform portion of just 4% (PDA) wh.ch agrees quite well with the lowest 

to Weld hi' Cn b> Stemer ’ Even thOU8h the 8auge threshold may not be expected 
to yield the same percentage in the model, th.s ind.cates the model hkelv 

underestimates the frequency of light rainrates as was found by Sui et al (1998) 

making it difficult to evaluate the performance of individua, methods blsld on 

rainfahthTt 1 7T Short et »'■ »»7) reported that 40% of the 

CCUrred dunng act,ve periods in TOGA COARE was stratiform in 
good agreement with the C & H, CCE and Xu methods in the model The In 

“ W 7 ld ^ h — A ” of *e separation techniques were dsi e 
m finding larger stratiform amounts in the tropical case. 


scheme irLlRKTORM aS caf U b d t t0 h" inVariant With res P ect fo advection 

EblORM case, but not so with the TOGA fOARP u 

— : ^ ll0W d Ce " S are pr ° duced leading edge of the squall system 

ng 4th order advection while PDA leads to fewer deeper cells that L , 

onger behind the leading edge resulting in more cells embedded with dn the 
::Zl re8,0n ' TWS " leadS * — between the ££ 


CA20 method was a reasonable approximation. “ ““ ” 


b o,h £%££?£££ coaIe“ ,he mean doud * 

». . — isr .tkt 

1 ' 0tal mass flu * Profiles for both types of advection were nearly identical a, 
all levels in the PRESTORM case bu, differed in the TOGA COARF ti 

shapes PRESTORM 

peak upward mass flux varied a great deal in the convective 



region (from 3.5 to 6.5 km for PDAi with u • • 

peak (between 8 and 10 km). I n the TOGA^COARE ^ * be anviI 

techniques showed much less • case ' the separation 

prohJal, peaked neTr 5 ^ "*“« ‘»e convective 

due to the lack of s.enifican, ve^r ^ ab °- * km 

edge. C&H GCE and CA9n • m ° tl0n ^ J ° W ,eVels behind the leading 

5 ^oin, and CA20 assumed nearlv all M g 

aloft to the stratiform region while increasingly Xu StetaerTnTvt w' 5 ^ 

them to the convective region significantly Altering the shape oTihe “T" 

Thus quahtatively the separation techniques were found o 1 ? 

significantly in terms of magnitude. q d ^ but differ 


advection fo “" d b e independent of 

result of the difference in C^str i "° l , S ° f ° r ^ C ° ARE A « ai " » a 
the PDA simulation while 4th order heat'"' 3 ' leat ' n8 W3S larger ab ° Ve 10 km in 
Agam, qualitatively the separati™ 5^^ *"* « ^ 
stratiform Ql pro fj] P c RoT^ a ^ ^ n terms convective and 

— i-*. Jz’Jszzsssrs ~:r r? 

stratiform reg.on while increasingly in Xu.^einm and WW T" * ^ 
to the convective reeion rp«.,lH«n. • Vt/W assi 8 n more 

peak anvil he.dpgTT.T^" “ ° Ve ”" ° f ° 2 *» ° « (deg/mm, of 

peak at a higher level and heating occurs ^ a'muchT’T'T hea ‘ inS Pr ° fi ' eS 
the tropical profiles. In the anvil, low-level coolinT d ^ C ° mpared t0 
pronounced in the tropical profiles while latent heating afofUs gTeal^ m ° re 

nearly t :r^^ is ,Q2) ^ runs are 

p ern is very similar with the TOP a 

except that convective drying is a little stronger for the 4th orf reSU ' tS 

hence the total profile as well Th . e 4th order simulation and 

C* Profiles is LeT, ^ ^ a " d 

COARE runs where the variation ,he , ( PRESTORM ««n in the TOGA 

largest convective drving and a I ^ ° ^ CA2 ° profile s have the 

account for the larger variation " ' m ° ,Stenin « in the ERESTORM case and 
ger tartation. The deeper cells in the TOGA COARE PDA run 
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do allow for some variation above 4 km where the Vt-W method shows the 
widest separation between convective drying and anvil moistening. Again all 
the methods show a qualitative consensus. 

Despite the fact the GCE method is the only one that directly uses 
hydrometeor information to differentiate between convective and stratiform 
areas, the amount of convective cloud water or cloud ice for the GCE method 
relative to the others followed the overall rank in convective-stratiform 
separation. All of the methods found the bulk of cloud water to be convective in 
the PRESTORM case. The results for cloud ice were split ranging from nearly 
equal parts convective and anvil with C&H to over 80% convective with Vt/YV. 
For TOGA COARE, the variation in the amount of convective cloud water was 
large while nearly all the methods expect for Steiner and Vt-W had cloud ice as 
overwhelmingly stratiform. 

Contoured Frequency with Altitude Diagrams (CFADs) were constructed 
for both the PRESTORM and TOGA COARE cases for each of the various 
separation algorithms. The PRESTORM reflectivity CFADs for each of the 
separation algorithms did indeed have larger peak modes in their convective 
CFADs. They also all had stratiform CFADs that showed evidence of a bright 
band. However, only Steiner and Vt-W had convective CFADs that were 
broader than their stratiform ones. The trends were similar for the TOGA 
COARE case except that the CA20 method produced somewhat broader 
convective CFADs and the stratiform CFADs for the PDA run did not show 
evidence of low-level evaporation. Vertical velocity CFADs for the PRESTORM 
case show that all of the methods have broader distributions in the convective 
region. The stratiform CFADs did show instances where points were clearly 
misclassified. CFADs of vertical velocities for the TOGA COARE case show 
similar behavior. Just as with PRESTORM, Vt-W produces the narrowest 
stratiform vertical velocity CFADs followed by Xu and Steiner. Thus all of the 
methods are able to satisfy the condition of having larger convective peak 

frequencies but not necessarily distributions that are broader for the convective 
regime. 


Modified CFADs were constructed on the microphysical processes in the 
model to show how the processes were distributed as a function of W and Z. 
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Most of the microphysical processes were distributed over a wider velocity and 
vertical range in the midlatitude case. The processes were also partitioned into 
convective and stratiform portions for each method. Condensation varied from 
98 (Steiner and Vt/W) to 71% (C&H) convective in the PRESTORM case and 92 
(Vt/W) to 55% (C&H) convective in TOGA COARE. Evaporation pretty evenly 
split between convective and stratiform in PRESTORM and was mainly 
stratiform in TOGA COARE. Deposition varied greatly ranging from 70 (Vt/W) 
to 31% (C&H) convective in PRESTORM to between 68 (Vt/W) and 4% (C&H) in 
TOGA COARE. Sublimation was clearly stratiform though Steiner and Vt/W 
divided it nearly evenly in TOGA COARE. Freezing was strongly convective in 
PRESTORM but varied considerably in TOGA COARE from 13% (C&H) to 92% 
(^t/W) convective. Melting mainly occurred in the convective region in 
PRESTORM while the opposite was true for TOGA COARE. The values are 
inclusive of both forms of advection. Thus depending on the process, the 
methods varied widely in their characterization and more so in the tropical case. 

Latent heating retrieval was shown to be quite sensative to the use of 
separation technique. This was mainly attributed to the stratiform region 
especially for methods that found very little stratiform rain (i.e., Steiner and 
Vt/W). 


Overall, it was found that the different separation techniques produced 
results that qualitatively agreed. However, the quantitative differences were 
significant. Observational comparisons were unable to conclusively evaluate the 
performance of the techniques. Thus it is important when using various 
convective-stratiform separation the sensitivities involved. Although not 
addressed in this study, it may be useful to add an additional category in the 
separation of convective systems: the transition region. Indeed, some of the 
differences between the methods studied here are related to their treatment of 
the transition area. C&H would be the most likely method to treat it as 
stratiform and Vt/W the most likely to treat it as convective. As cells move 
from the leading edge of the system rearward, it can be very arbitrary when to 
change their designation from convective to stratiform. As such, based on the 
fact that many of the separation techniques look for convective criteria only 
below the melting level, it may be useful in the future to define the convective 
region as that where convective criteria are exceeded below the melting level, the 



transition region as that where convective criteria are only exceeded above the 
melting level, and the remaining region as stratiform. 
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TABLES 


Table 1 Characteristics of the large scale environments associated with the 
PRESTORM and TOGA COARE squall line cases. Simulation domain 
size and grid spacing and references are also listed. 


Table 2 Main criteria used by the six different separation algorithms to partition 
convective and stratiform regions. 


Table 3 


Total accumulated rainfall (mm) per grid point over all 12 hours of 
simulation time for all four PRESTORM and TOGA COARE 
simulations. Mean stratiform percentages for all six convective- 
stratiform separation techniques are also shown for each simulation. 
Values in parentheses are for the last 4 hours of simulation time. 


Table 4 Mean width (km) of the convective and stratiform regions over 12 hours 

of simulation for all six separation methods and their average for all 
four simulations. 


Table 5 Total mean percentage of the microphysical processes that are convective 
according to each separation technique for the PRESTORM simulation 
with 4th order advection. 

Table 6 Total mean percentage of the microphysical processes that are convective 

according to each separation technique for the PRESTORM simulation 
with PDA. 


Table 7 Total mean percentage of the microphysical processes that are convective 
according to each separation technique for the TOGA COARE simulation 
with 4th order advection. 


Table 8 Total mean percentage of the microphysical processes that are convective 

according to each separation technique for the TOGA COARE simulation 
with PDA. 



FIGURE CAPTIONS 


4 7 


Fig. 1. Vertical cross-sections of model estimated radar reflectivity (dBZ) after 
720 minutes of simulation for (a) the PRESTORM case with 4th order 
advection, (b) the PRESTORM case with PDA, (c) the TOGA COARE case 
with 4th order advection, and (d) the TOGA COARE case with PDA. The 
convective regions for each of the separation techniques are overlayed in 
solid black lines. Traces of the corresponding surface rainrates (mm/h) 
are also overlayed (top of the domain = 200 mm/h). 

Fig. 2. Time-domain cross-sections of surface rainfall rates (mm/h) for (a) the 
PRESTORM case with 4th order advection, (b) the PRESTORM case with 
PDA, (c) the TOGA COARE case with 4th order advection, and (d) the 
TOGA COARE case with PDA. 

Fig. 3. Time series of instantaneous grid averaged total (solid) and stratiform 
(dashed) rainrate for both types of advection for the (a) PRESTORM and 
(b) TOGA COARE cases. 

Fig. 4. Vertical profiles of mean cloud coverage as a percentage of the 
horizontal domain for both types of advection for both the PRESTORM 
and TOGA COARE cases. 

Fig. 5. Vertical profiles of domain average accumulated mass flux for the total 
(heavy dash), convective (solid) and stratiform (dotted) regions for both 
types of advection for the (a) PRESTORM and (b) TOGA COARE cases. 

Fig. 6. Vertical profiles of domain average accumulated convective and 
stratiform mass flux for each of the six separation techniques using PDA 
for the (a) PRESTORM and (b) TOGA COARE cases. [C&H - black dashed, 

GCE - solid black, Xu - gray dotted, CA20 - solid gray, Steiner - black 
dotted, Vt-W - "#"] 

Fig. 7. Vertical profiles of rain normalized domain average apparent heating 
(Ql) for the total (heavy dash), convective (solid) and stratiform (dotted) 
regions for both types of advection for the (a) PRESTORM and (b) TOGA 
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COARE cases. 

Fig. 8. Vertical profiles of rain normalized Ql components for the (a) 
PRESTORM convective region, (b) PRESTORM anvil region, (c) TOGA 
COARE convective region, and (d) TOGA COARE anvil region for PDA. 
[Condensation - solid. Deposition - thin solid. Evaporation - dotted. 
Sublimation - thin dotted, Freezing/Melting - thick dashed, Eddy flux - 
"e ", Radiation - "+"] 

Fig. 9. Vertical profiles of rain normalized domain average convective and 
stratiform Qi for each of the six separation techniques for (a) the 
PRESTORM case with 4th order advection, (b) the PRESTORM case with 
PDA, (c) the TOGA COARE case with 4th order advection, and (d) the 
TOGA COARE case with PDA. Line patterns follow Fig. 6. 

Fig. 10. Vertical profiles of rain normalized domain average apparent 
moistening (Q 2 ) for the total (heavy dash), convective (solid) and 
stratiform (dotted) regions for both types of advection for the (a) 
PRESTORM and (b) TOGA COARE cases. 


Fig. 11. Vertical profiles of rain normalized domain average convective and 
stratiform Q 2 for each of the six separation techniques for (a) the 
PRESTORM case with 4th order advection, (b) the PRESTORM case with 
PDA, (c) the TOGA COARE case with 4th order advection, and (d) the 
TOGA COARE case with PDA. Line patterns follow Fig. 6. 

Fig. 12. Vertical profiles of domain average hydrometeor content for (a) the 
PRESTORM case with 4th order advection, (b) the PRESTORM case with 
PDA, (c) the TOGA COARE case with 4th order advection, and (d) the 
TOGA COARE case with PDA. [Rain - solid black. Cloud - dotted black, 
Snow - solid gray, Cloud ice - dotted gray. Hail - thick dashed black, 
Graupel - thick dashed gray] 

Fig. 13. Reflectivity CFADs for the (a) PRESTORM convective region using GCE 
separation, (b) PRESTORM anvil region using GCE separation, (c) 
PRESTORM convective region using Steiner separation, (d) PRESTORM 



anvil region using Steiner separation, (e) TOGA COARE convective 
region using GCE separation, (f) TOGA COARE anvil region using GCE 
separation, (g) TOGA COARE convective region using Steiner 

separation, (h) TOGA COARE anvil region using Steiner separation for 
PDA. 


Fig. 14. Same as Fig. 13 except for vertical velocity and Vt-W separation in place 
of Steiner. 

Fig. 15. Domain total microphysical CDADs of (a) PRESTORM condensation, (b) 
PRESTORM evaporation, (c) PRESTORM deposition, (d) PRESTORM 
sublimation, (e) PRESTORM freezing, (f) PRESTORM melting, (g) TOGA 
COARE condensation, (h) TOGA COARE evaporation, (i) TOGA COARE 
deposition, ()) TOGA COARE sublimation, (k) TOGA COARE freezing, (1) 
TOGA COARE melting, 

Fig. 16. Simulated profiles of retrieved heating (Qi) for each of the six separation 
techniques applied to model data (PDA) but using observed stratiform 
percentages for (a) PRESTORM and (b) TOGA COARE. The actual model 
heating is shown by the thick solid black line. Line patterns follow Fig. 6. 

Fig. 17. Convective and stratiform heating profiles normalized by their 
respective rainfall amounts for each of the six separation techniques for 
the (a) PRESTORM convective region, (b) PRESTORM anvil region, (c) 
TOGA COARE convective region, and (d) TOGA COARE anvil region. 



Table 













Table 





Table 




Table 






Table 


X 

4-> 

s 

PH 

o 

H 

C/5 

W 

P£5 

Pin 


ID 


Vt-W 1 

97.6 

76.5 

C 

IN 

41.5 

96.7 

81.3 j 

Sh 

a> 

C 







CN 

CO 

CO 

GO 

00 

<N 

*3 

«w 

IN 

c\i 

rsi 

CO 

CO 

d 

ON 

tN 

CO 

CO 

ON 

00 

CO 







o 

<N 

IN 

IN 

rH 

<N 

in 

in 

< 

GO 

IN 

rH 

in 

as 

rH 

u 

GO 


rt< 

<N 

tN 

in 

3 

in 

ON 

rH 

IN 

rH 

ON 

x 

c4 

in 

CO 

in 

n! 

rH 


Os 

m 

in 

<N 

GO 

CO 

ac 

CO 

CO 

IN 

m 

O 

in 


d 

CO 

rH 

• 

rH 

in 

04 

u 

IN 


CO 

<N 

co 

rH 

w 

CD 

rH 

CO 

rH 

GO 

00 

u 

CO 

00 

CD 

04 

o 

rH 

o 

00 

rf 


04 

00 

in 


c 

0 

1 pH 

«4-» • 

c 

0 

1 

c 

o , 

c 

( 0 

to 

60 


C3 ’ 

C/5 

c 

OJ 

-u> • 

c3 ; 

Vh 

o 

s 

f— 1 

Cft 

o , 

£ 

c 

N * 
0> * 

C 

-M 

* 

xt 

Dh 

Cu ’ 

1 ^ 

<u , 

0/ 

( 

c 

5 1 

c3 

> 

Jj 

0> , 

^ « 

3 [ 
J) 

u % 
in * 

s 




Table 6 PRESTORM/PDA 







< 

Q 

Ph 

w 

& 

< 

O 

u 

< 

o 

o 

H 

oo 


H 


















































PRESTORM 



(g/kg/hr) / (mm/hr) 


TOGA COARE 



(gAg/hr) / (mm/hr) 










HEIGI- 









XJ PRESTORM pda 

CFAD - Refl 
U 1 1 1 1 1 1 1 1 1 1 1 i 1 m i i i ) i i r i i i r t 
I Stratiform - GCE 



(dBZ) 


PRESTORM pda 

CFAD - Refl 

L T i i f i ri - ! ! 1 1 i i i i i i i | | | | | | , r ri ~j 

- Stratiform - Stn r 






TOGA COARE pd 


100L Convective 



TOGA COARE pda 
CFAD - Refl 


100L Stratiform — 




g 

~i' i i i f r 

TOGA COARE pda 
CFAD - Refl 


h 

TOGA COARE pd 
CFAD - Refl 

Lvective 

-- Stnr 

1 « ' i 19 

■ 

1 I 1 T“7 1 r*T — ! 

100 ~ Stratiform 

r 1 1 1 1 1 1 1 1 1 1 ITTT 

1 -1C 

— Stnr 





ffl 


E 









































100L TOTAL 


TOGA COARE pda 
CDAD - Cond(W) 



TOGA COARE pda 
CDAD - Evap(W) 


loot TOTAL 


0 5 10 15 



loot TOTAL 


TOGA COARE pda 
CDAD — Depo(W) 



loot total 


TOGA COARE pd 
CDAD - Subl(W) 



SFC UJ - J - J 1 
-10 -5 


TOGA COARE pda 
CDAD - Frez(W) 



rOGA COARE pda 

CDAD - Melt(W) 














HEIGHT(km) 











